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Abstract 



We study the entanglement entropy of a scalar field in 2 + 1 spacetime where space 
is modeled by a fuzzy sphere and a fuzzy disc. In both models we evaluate numerically 
■ the resulting entropies and find that they are proportional to the number of boundary 

degrees of freedom. In the Moyal plane limit of the fuzzy disc the entanglement entropy 
per unite area (length) diverges if the ignored region is of infinite size. The divergence is 
(interpreted) of IR-UV mixing origin. In general we expect the entanglement entropy per 
unite area to be finite on a non-commutative space if the ignored region is of finite size. 

in 
o . 

1 Introduction 

It is believed that black hole thermodynamics hold important clues to the nature of the 
structure of spacetime at the very small scale . For example the finiteness of black hole entropy 
is understood as a direct manifestation of the discreetness of spacetime at the Planck scale, and 
points out to a necessary reduction of the number of degrees of freedom on the horizon [1,2]. 
A piece of evidence for this comes from the contribution of the thermal atmosphere ( i.e the 
entropy stored in the field quanta near the horizon ) to the black hole entropy which was 
shown to be divergent in absence of any UV cutoff which goes against the finiteness of black 
hole entropy [3]. This divergence is closely related to the entanglement entropy ( or geometric 
entropy ) considered in [4,5] and later in [6]. 

The notion of entanglement entropy originates from the simple observation that an observer 
outside the horizon has no access to the degrees of freedom behind the horizon. For this 
reason the outside observer would describe the world with a reduced density matrix obtained 
by tracing out the unaccessible degrees of freedom behind the horizon. If the exterior modes 
and the external modes are correlated "entangled" the resulting density operator is thermal 
even if the global state of the system were pure. 

The entanglement entropy has a UV divergence which can not be renormalized away, at 
least when calculated for fields in a fixed background. Moreover in [8] it was shown that the 
divergences in the entanglement entropy are the same divergences one must deal with when 
trying to renormalize the theory of gravity coupled to matter, and therefore understanding the 
finiteness of the entropy may not be possible without a complete knowledge of the UV behavior 
of the theory. 
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However if we introduce a short distance cutoff around the planck length one obtains a finite 
entropy of the same order of magnitude as that of a black hole. This entanglement entropy 
generally scales like the area of the boundary. This has led to many speculations attributing 
the black hole entropy to the sum of all entanglement entropies of the fields in nature [11]. 
Whether or not the entanglement of quantum fields furnishes all of the entropy or part of it, 
contribution of this type must be present, and any consistent theory must provide for them in 
its thermodynamic accounting. 

Although the entanglement entropy has been extensively explored by many authors most of 
the work was oriented toward computing it using different regularization methods and studying 
the degree and the universality of the divergences [7,9,10]. For instance in [9] it was shown 
that except in 1 + 1-dimensions the coefficient of the divergences is non-universal and depends 
on the regularization method. Therefore any quantitative comparison of different calculations 
was not possible. 

The aim of this paper is to study the entanglement entropy on some fuzzy spaces and draw 
possible conclusions about the entanglement entropy on non-commutative spaces in general. 
The motivations for considering the entanglement entropy on fuzzy spaces are many. One 
obvious reason is the general hope that non-commutativity would soften UV divergences and 
may render UV divergent quantities finite as it is the case on fuzzy spaces, and at the same 
time preserve the basic symmetry of the space . For instance results from finite-temperature 
non- commutative field theory showed that the non- commutative model behaves as if it had 
many fewer degrees of freedom in the UV than in conventional field theory. In some cases 
the degrees of freedom were so drastically reduced that the UV catastrophe could be avoided. 
See [12] and references therein . 

Another motivation is that non-commutative geometry arises naturally in some limits of 
string theory in connection with D-branes. In effect it has recently been suggested by many au- 
thors that the microstates counting of non extremal black holes using field theory dual of string 
theory could be interpreted as arising from entanglement [15-17]. In considering entanglement 
entropy on fuzzy or noncommutative spaces one faces the question of defining the boundary 
with respect to which the entanglement is to be computed. Indeed in view of the non locality 
of non commutative theory one would expect any boundary to be fuzzy. This also brings in the 
question how to define disjoint regions in non-commutative or fuzzy spaces in order to properly 
define the entanglement entropy. We will show in this paper that once the field variables are 
properly chosen there is almost a unique way to obtain entanglement between different degrees 
of freedom in different regions. 

This paper is organized as follows. In section 2 we give a brief review of the formalism 
which we will use. In particular we compute the entanglement entropy for a free scalar theory 
defined on the continuum sphere. In section 3, part 1 we compute the entanglement entropy 
resulting from tracing half of a fuzzy sphere and show that the result is proportional to the 
number of boundary degrees of freedom or equivalently to the area of the boundary . In part 
2 of section 3 we reconsider the same problem on a truncated Moyal plane, i.e a fuzzy disc, 
and obtain similar results as in fuzzy sphere case. We also discuss the Moyal plane limit and 
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observe that if the ignored region is blown up a UV-IR mixing phenomenon takes place. We 
conclude in section 4 by general discussions and possible improvement of our numerical results. 



2 The Regularized theory 

In this section we review the main formalism we will be using and apply it to a regularized 
continuum sphere. For detail of the formalism we refer the reader to [5,6]. For a review see [13]. 
The Hamiltonians that we will consider in this paper are of the standard form 

H=\ Y,^a,bit a k b + VW V) (2-1) 

A,B 

Vab is a real symmetric matrix with positive definite eigenvalues. The case where V has 
zero eigenvalues needs special treatment. It corresponds to the case of a massless field and 
we will return to it later. The normalized ground state of (j2.1|) is given in the Schrodinger 
representation by 

< ip A | >= [det — ] 1/4 exp [ - ^WWV] (2-2) 
Where W is the square root of the matrix V. The corresponding density matrix is 

p{<p,<p') = [det-] 1/2 exp[- V + ^'V*)] (2.3) 

7T Z 



Now if we consider the information on the fields degrees of freedom <p a , a = 1, n as unavailable, 
we form a reduced density matrix by integrating over the (p a , a = 1, n, 

/n 
n^v^v) (2.4) 
a=l 



The entanglement entropy is the associated Von Newman entropy of p re d defined by S = 
— Trp re( i logp re d. The entanglement entropy for any Hamiltonian of the form ()2.1|) can be 
shown to be given by 



(2.5) 



Sent = l0 § + Vl + Ailog + ^1 + ^ 

Where Aj are the eigenvalues of the following matrix 

n 



a=l 



W a j and W^ 1 are elements of W and W 1 respectively with i, j running from n + 1 to N and 
a from 1 to n . A is an (N — n) x (N — n) matrix and i, j run now from n + 1 to N. 
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Before considering theories denned on fuzzy spaces we start with a real free scalar field 
defined on spacetime with topology R x S 2 where R is time and S 2 the spatial slice. The 
lagrangian is given by 



L = l - I dn( 12 



(A-^ 2 )0). 



(2.7) 



For computational reasons it turns out that the theory is easily regularized in the cylindrical 
coordinates. The Laplacian A in this coordinates is given by 

Id d Id 2 

A = R?!h^ R ~ Z W + R 2 - z 2 d(t> 2 ' 

Expanding the field in Fourier modes by writing 



(2.8) 



j +00 

;== Y] <Pm{.z)e 
Ztt 



-irmp 1* 
5 Ym 



Z = 0-m.{Z) 



(2.9) 



We find that the Lagrangian L is the direct sum of individual microscopic Lagrangians L m each 
associated with an allowed value of the azimuthal number m, viz 



where 



R 



dz 



Q 



L= J2 L r. 



1 _ — \( \ 2 _ i 

B 2>[ ~dz~ ) ~ [ 



(2.10) 



m 



R 2 -z< 



+ V 2 )Q 2 m 



(2.11) 



where Qq = 4>(z), Q m = \f2Re(j) m for m > and Q m = y/2Im(J) m for m < . 

We regularize this model as follows. The z— axis is replaced by a one- dimensional lattice, 
i.e z — >z n = na where a = ^ — >0 is the lattice spacing and n = — (iV — 1), .., [N — 1). The 
regularized microscopic Lagarangian for one sector is given by 



1 N ^ 
Lm = Ya 22 

n=-(7V-l) 

The corresponding Hamiltonian is 



Th 2 Tft 

aQ m , n ) 2 - (1 - jp) (Qm,n Qm,n—l) ^j\T2 _ n 2 ^ ^Qm, 



(2.12) 



N-l 



1 E 

2a ^ 



7T7 



Th~ 2 
t n + (1 — ) [Qm,n ~ Qm,n—l) 



m 



N 2 -n 2 



(2.13) 



n=-(JV-l) 

By scaling the fields and shifting the summation variable the Hamiltonian H m can be brought 
to the general form of equation (j2.1|) . namely 

2N-1 r 



E 

A,B=1 



6 A>B 7T A ^ + VXbQzQ 
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(2.14) 



where 



V ab - (2 — — +______ 

(N-A)\ x (N — B) 2 
- (1 ^2 )<W-i-(l ^ (2.15) 

The macroscopic Hamiltonian is given by 



+oo 



H= H™ (2.16) 



m=— oo 



If we now consider the degrees of freedom residing in the upper hemisphere unaccessible we 
construct the reduced density operator for the ground state by integrating all the modes 
for a = N, 2N — 1 corresponding to positive z for all values of m. According to (J2.16j) the 
resulting reduced density operator is of the form 

oo 
m=— oo 

with pj™^ = Pr ed m \ total entanglement entropy is therefore given by 

oo 

S N = S + 2j2S m (2.17) 

771=1 

S m is computed by applying equation (|2.5j) to the Hamiltonian H m with 

2N-1 



a=N+l 

The indices i and j run in the available region i.e i,j = 1, ...,N. For a given iV we compute 
numerically the entropy for different masses. This entails the computation of the eigenvalues of 
A for each pair of values (N, m) and then computing the entropy by means of equation (J2.17I) . 
However to obtain a numerical result the sum over m must be cutoff at some value m max . This 
can be decided by the following simple observation similar to that in [6] . For large m, m 3> N, 
the m-dependent term dominates over the other terms and S m can be computed perturbatively. 
It is found that for m^> N 

A^ = ^^^ + 0((-) 6 ) 
4m 4 m 

which gives for S m 

This demonstrates that the sum over m will converge and also can be used to set an upper 
bound on the remaining of the sum which can easily be seen to be negligible for m ^> A 4//3 . 
The large m behavior of S m is very similar to that obtained in 3 + 1 by Srednicki [6]. 
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Figure 1: The scaled entropy S for a massless scalar field on continuum sphere, with m r 
2000. 



The result for the scaled entropy S = S^/N, m max = 2000 and massless field is depicted 
on figure 1. It is observed that the entropy converges for large N to the value 0.465iV. Hence 
the resulting entropy is proportional to the area (circumference) of the boundary with a pro- 
portionality coefficient which goes like 1/cutoff, via 

S, = °-^A (2.20) 
where A = 2tcR, which is the area law in 2 + 1 dimension. 



3 Fuzzy spaces 

3.1 The fuzzy sphere 

We start with a scalar field on RxS^, where S 2 tv is a fuzzy sphere of matrix dimension 
N = 21 + 1. The action reads 



S N = ^JdtL, L^^TrU 2 + <f>[C 2 -^]A (3.1) 

The scalar field is an NxN hermitian matrix with mass parameter /i. The Laplacian £f is the 
SU (2) Casimir operator given by C 2 = C\ + C\ + £| with action defined by Ci(<f>) = [Li, 0] and 
C 2 ((f>) = [Li, [Li, <j>]\. The Li satisfy [Li,Lj] = ie^Lk and they generate the SU{2) irreducible 
representation of spin I = ^f^. 
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First we observe that 



,cd4>cd (3.2) 

where 

M ah)Cd = ^{Ll) hc 5da — {.L,^) bc (Li) d ^j = ^c 2 5b c Sda — (Li) bc (Li) d ^j . (3.3) 

c 2 is the Casimir of the spin / IRR of SU(2), i.e c 2 = 1(1 + 1). Now we introduce two real 
matrix scalar fields by splitting as follows = <f>i + 20 2 . The hermiticity of implies that 
<Pi = 0i and <p2 = —02, in other words 0i is a real symmetric NxN matrix whereas 02 is 
a real antisymmetric NxN matrix. By using also the fact that the matrix M abtCd satisfies 
Mba,dc = M abjCd it is a straightforward calculation to show that 

Tr0£^0 = 2(0x + 2 UM aMc (0x + 2 ) cd = 2$ ab M aMc $ cd . (3.4) 

where $ = 0i + 2 . For explicit calculation we use the matrix form of L\ and L 2 in the base 
where L 3 is diagonal. They are given for arbitrary I by 

1 % 
{L\)ab = -[Bb5 a ,b+i + B a 5 ajb -i] , (L 2 )ab = - [Bb5 a ,b+i — B a 5 a ^-i] , (L 3 ) ab = A a 5 a:b (3.5) 

Where B a = ^a(N - a) and A a = -a + The indices a, 6 run from 1 to N = 21 + 1. We 
can immediately compute 

Tr0£^0 = 2 $ a5 (c 2 - A a A 6 )$ a6 - ^(So-iSft-i)^.!,^! - <£> ab (B a B b )$ a+hb+1 . (3.6) 

Now, this expression suggests the introduction of the following fields : defined by Q a m ^ = 
&a,a+m and Qa~ m ^ = & a +m,m for m = 0, ..N — 1. More explicitly we have for m = the field 
Q^ -* with N degrees of freedom given by 

Q(0) = ($11,^2...^ (3J) 

For m = +1 we have the fields and each with N — 1 degrees of freedom given by 

Q(+l) = ($12 ) $ 2 3) $34) ^ $N-l,Nj ? g(-l) = $ 32 ? $ 43 ? „ ^JV-l) (gg) 

For general positive m we have the fields Q < - +m ' ) and Q^'^ which contain each (N — m) degrees 
of freedom given by 

The last two fields and Q - ^ -1 ) contain one degrees of freedom each, viz 

Q(N-l) = ($ liY } ? Q (-N+l) = (QNly ( g 1Q) 
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Using this parametrization we can show that 



JV-l N-\m\ 

TrjCfy = 2 E <?< m) 

m=-(N-l) a=l 



(C2 — ^4a^a+|m|)^a,6 ~ 2^o-l-Bo-l+|m|^o-l,6 — ^ ^ a ^ a +l m l^ a+1 ' 6 



Similarly we can compute 



7>(o 2 - /r 



JV-l JV-|m| 

)= E E Q^i-dl + ^Q^- 

m=-{N-\) o=l 



(3.12) 



Hence the Hamiltonian if of the free theory takes the form 

JV-l N-l N-\m\ 



h= E H m= E E 



m=-(JV-l) 



m=-(JV-l) a,6=l 



(3.13) 



where 



v (m) = 2 



(C2 + "2 A a A a+ \ m \)5 a £ — --B a -l-Ba-l+|m|^a-l,6 — ^ B a B a +\m\ S a +l,b 



(3.14) 



(m) 



Ofa™^ 1 ■ With this result one can see that the free theory splits into 2(2/) 



and ii. 

independent sectors {7i m },m 



-(N 



•, (N — 1), each sector 7i m has N — \m\ degrees 



of freedom (N — \m\ coupled harmonic oscillator) and described by a Hamiltonian H„ 
ground state density matrix is easily seen to be 



JV-l 
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(m) 



The 



(3.15) 



-(N-l) 



It is worth mentioning here that the formalism of the fuzzy sphere bears a lot of similarity to 
the lattice regularization of the continuum discussed in the previous section. However the fuzzy 
sphere provides a natural cutoff for the quantum number m which is playing a similar role to 
the azimuthal number of the continuum sphere. 

Having brought the Hamiltonian on the fuzzy sphere to the form we want we are ready now 
to discuss the entanglement entropy in the fuzzy sphere setting. First we note that in order to 
introduce the entanglement entropy we need to divide the field degrees of freedom into two sets 
residing in two disjoint regions corresponding, say, to the upper and lower fuzzy hemispheres. 
In order to do that rigourously one needs to give precise criterion that allows one to decide 
whether a given set of two fields ( two matrices) have a disjoint support on the fuzzy sphere. 
Although we believe that deriving such criterion may not be difficult we shall content ourselves 
to some heuristic and intuitive arguments. We first note that the form of the Hamiltonian 
suggests strongly that we take each sector H m and trace over half of the degrees of freedom. 



Q { b m) - 

(3.11) 
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For a fixed N and m the number of degrees of freedom in the sector 7i m is N — \m\, UN — \m\ 
is even we trace out the following degrees of freedom 

Q { r\Qt\---iQ { r\k = ^^ ( 3.i6) 

if N — \m\ is odd we have two options, either we trace out 

Q?\Q? , ,--;Qt ) ,k= N -' i f- 1 (3.17) 

or we trace out 

QS m, .gr ) ,---,,gr,t-= jV ' l 2 m|+1 (3.18) 

However both options lead to the same entanglement entropy for large N and the degrees of 
freedom Q ™-\ m \+i wm be interpreted as boundary degrees of freedom and there are N of them. 

2 

This corresponds in the original matrix notation to dividing the matrix into two parts, upper 
left triangle (f>u and right lower triangle <pL ■ For example for the first option above for N = 5 
the <f)jj and 4>l will look as follows 



Hj 



/ 


011 


012 


013 


014 


o\ 




/ 














015 


\ 




021 


022 


023 





















024 


025 






031 


032 











, 0L = 










033 


034 


035 






041 





















042 


043 


044 


045 




V 














°/ 




V 


051 


052 


053 


054 


055 


/ 



(3.19) 



The components 05i, 042, 033, 024, 015 are the boundary degrees of freedom. (f)jj and 0£ can 
be given the interpretation of corresponding to functions with disjoint supports, one on the 
lower half and the other on the upper half of the fuzzy sphere. A hint for this comes from 
the observation that Tr(0[/0 £ ) = 0, for arbitrary 0^ and 0^ and the fact that we are working 
in the basis where L 3 is diagonal and therefore one can talk about negative and positive z 
coordinates ( i.e upper and lower hemisphere). Indeed this is the only choice that leads to a 
non-zero entanglement entropy. Now from equation (J3.13j) one can easily see that the resulting 
reduced density operator will be 



N-l 



Pied 



(m) 
Pred 



(3.20) 



m= _(jV-l) 



and the associated entropy is 



21 



3 l = 3 + 2^2s m , N = 2l + 1. 



(3.21) 



m=l 
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where S m is computed using equation (j53J) with A< m ) defined by A[f = - £* =1 w ^W aj where 
W is the square root of the potential matrix given by equation (j3.14j) . The indices i,j 
run from k + 1 to iV — m (the available region). Since we expect our entropy to be proportional 
to the circumference of the equator which can be said to be a fuzzy circle, we define a scaled 
entropy by dividing the entropy by the square root of the Casimir C2, viz 

S = - Sn (3.22) 

We compute the scaled entropy S numerically for N = 4 , 200 and run several numerical 

calculations with different values of the mass. The results are depicted on figure 2 . We can 
immediately see that the different curves with masses /x 2 = 0, 10 -5 , 1~ 3 , 1 approach the same 
value for large N. For example for I = 600 the scaled entropy for the different masses is equal 
to 0.39 and the discrepancy between the four different values being in the third digit. 

We can thus safely extrapolate to very large / or N and conclude that the entropy is given 

by 

S N = 0.39^1(1 + 1). (3.23) 

Since the radius of the fuzzy sphere can be scaled with I keeping the noncommutativity param- 
eter 9 = , R fixed, we can write the entropy in the following form 

y/ 1(1 + 1) 

Sl = 0.39? = ^ (3.24) 

This is exactly the area law in two dimension. 

Some remarks about our numerical evaluations for the entropy are in order. As in any 
numerical calculation, numerical errors are inevitable. In our calculation numerical errors orig- 
inate from rounding errors and iteration, however it turns out that they are really tiny, in 
worst cases of order 10~ 10 . Indeed none of our matrices is ill-conditioned or nearly singular to 
a working precision, therefore one should not for instance expect significant errors to occur by 
inversion. 

Now, despite the fact that numerical errors are really small the exact numerical value for the 
leading contribution to the scaled entropy can not be extrapolated from our numerical result, 
because as figure 2 shows there are subleading corrections which are nonvanishing, although 
small, for finite I and affect the rate of convergence. 

Although one can not obtain the form of the subleading corrections, our numerical calcula- 
tions show that they are negligible for large I, for example they are of order 10 -3 for I = 600 
and for the values of the mass we used. 

Having established the area law, let us now turn to a simpler and more interesting inter- 
pretation of the result obtained. Since equation ()3.23|) is valid for large / we can write it as 

Si = 0.19(2/ + 1) ~ log2 2m (3.25) 

Therefore the entropy is directly proportional to the number of the degrees of freedom on the 
dividing boundary, in this case being a "fuzzy circle". This result fits with the picture of the 
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Figure 2: The scaled entropy 5* resulting from tracing out half of the fuzzy sphere. 

horizon (simulated here by the equator) as being divided or quantized into small cells of Planck 
size with each cell carrying roughly one bit of information [1]. The picture geos along the 
same line with old suggestion that the area of a black hole have a quantized spectrum. See for 
example [18]. 

We conclude this section by few comments about the result for the massless case. For 
relatively small values of N or / there is a clear mass dependence which dominates the behavior 
of the entropy. The entropy is a decreasing function of the mass as can be seen from figure 
2. However it shows a discontinuity at zero mass where it drops from large values for small 
masses to small values for zero mass. However this discontinuous behavior disappears for very 
large I and the entropy has a smooth dependence on the mass. To understand the zero mass 
discontinuity we go back to the original potential. It is easy to see from equation (|3.14|) that the 
potential of the sector TCq is a singular matrix in the massless case and therefore one can 
not apply directly the formalism of section 2. However, once treated carefully one can easily 
show that the sector 7i gives zero contribution to the entanglement entropy and therefore must 
be projected out from the sum in equation (|3.21|) . In the massless case the associated entropy 
is therefore 

N-l 

S N = 2J2 S m (3.26) 

m=l 

All the terms S m are computed using the same formalism of section 2 since none of the potentials 
V m is singular. Moreover, for large / (the continuum limit) the main contribution to S comes 



0.8 



<s> 0.6 



0.4 
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from the sectors with m different from zero, and the term Sq becomes irrelevant which explains 
the smooth behavior for large I. For instance already at / = 400 and for /x 2 = 0.1 both equations 
(I3.26j) and (|3.21|) give for Sn the same value S = 0.39 . Now the reason for which previous 
lattice computation, like the one presented in the first section and the one in [6], were not 
sensitive to this effect is due to the fact that those lattice regularizations break the translation 
symmetry or rotation symmetry to which the zero modes of the potential are related. 



3.2 Moyal plane and the fuzzy disc 

We consider now a scalar theory on RxRg where Rf> is now the Moyal plane. The action is 
given by : 

S = \J dtTr(<p 2 + 0(V 2 - /i 2 )0). (3.27) 

The trace is infinite dimensional and the Laplacian is given in terms of creation and annihilation 
operators a and a + by the experssion 

V 2 := -^[a + , [a, 0]] = -l[a, [a+ 0]]. (3.28) 

Let us recall that [a, a + ] = 9 where 9 is the noncommutativity parameter. The fuzzy disc is 
obtained from the plane following [14] as follows. We consider finite dimensional NxN matrices 
<f), viz 

N-l 

= ^ mn I m >< H I ' ^ + = ^ ' Cm = <t>mn- (3-29) 

n,m=0 

Then it can be shown that the Laplacian V 2 acts on a finite dimensional space of dimension 
(iV+ l) 2 , i.e V 2 is an(iV + 1) x (N + 1) matrix. The action on RxD^ is thus given by (l3~27l) 
where the trace Tr is simply cut-off at N. We denote this trace by Tr jy. The radius of the disc 
is given by 

R 2 = N9. (3.30) 

In other words when iV — >oo we must take 9 — >0 to recover a commutative disc. 

For the purpose of computing the entanglement entropy on this space the relevant piece of 
the action is given by 



4 8 8 

V = -Tr N (j)(V 2 -^)(j)=(iJ 2 + -)Tr N (f) 2 + —Tr N (j) 2 a + a-—Tr N (f)a(j)a + . (3.31) 



9 2 



9 2 



Explicitly we can compute 



iV-l 



n,m=0 



fi 2 9 



(1 + — + n + m)(f) nm (j) mn - v/^Ww^m-in-i + y/(n + l)(m + l)0„ m m+ i„+i 

(3.32) 
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or equivalently ( with <p nrn = (f) n -i m -i ) 



N 



n,m=l 



-1 + — + n + m)(f) nm (f) mn - y/(n-l)(m- l)<f> n m<f> m— In— 1 



(3.33) 



By using the same trick we used on the fuzzy sphere; namely we split the fields into symmetric 
part 4>mn = <j>nm and antisymmetric part 4>mn = —<f>nm by writing (4>) mn = 4>mn + i(/>mL, then we 
recombine into the real field <3> mn = 4>mn + 4>rnn, we can put the above action into the form 



e 



N 



n,m=l 
N 



- -Tr N 0(V 2 - t?)4> = 2 J2 *»».(- 1 + -J- + A « + A 'M 



E 

n,m=l 



n-1 n-lm-1 n+lrrt+l 



(3.34) 



where A n and B n are now defined by A n = n and B n = y2n. The off diagonal elements 
have exactly the same structure as the off diagonal elements on the fuzzy sphere, whereas the 
diagonal elements here involve the sum A n + A m as opposed to the product —A n A m on S^. 
Following the same steps we have taken in the fuzzy sphere case we can write the above action 
in the form 



-^Tr^(V 2 -^)0 



N-l N-\m\ 

E E vMqmqP 

m=-(N-l) a,b=l 



(3.35) 



where 



(m) 



ab 



2 
2 



fl 2 9 2 1 1 

( - H — ^ — \~ A a + A a+ \ m \)5 a)b - -S a _ 1 S a _ 1+ | m |5 a _i ;6 - -B a B a+ \ m \5 a+1;b 
(2a + \m\ - 1 + — )8 a ,b ~ a/(o - l)(a - 1 + \m\)5 a -i,b ~ V a ( a + H)<Wi,& 



(3.36) 



The fields Qj™ and Qi" m) ( m>0 ) are defined in the same way as on S^, in other words 



2Q^ = 2$ a , a+m = (1 - i)&-l,«*m-l + (1 + i)C-l, ffi+m -l 
2 Q(-m) = 2 $ a+mia = (1 - ^)0 a+m _ 1 , a _ 1 + (1 + i)^^. 



(3.37) 



Altogether we obtain 2iV — 1 independent sectors H m ,m = — N + 1, — iV + 2, TV — 2, N — 1. 

For the fuzzy disc we shall consider two different entanglement entropies. The first one will 
result from tracing out a smaller sub-disc, and the second one from tracing out half of the 
fuzzy disc. In the first case the ignored region will remain finite once we consider the Moyal 
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plane limit, while in the second case the ignored region blows up. Consider now the following 
(n + 1) x (n + 1) sub-matrix with n < N 



/ 


000 


001 


002 •• 


0On \ 




010 


011 


012 •• 






020 


021 






V 


0nO 






0nn / 



The degrees of freedom in this matrix have support on a sub-disc D n . We are going to assume 
that this region is unaccessible to outside observer and therefore we are going to trace the 
degrees of freedom residing in this region. One can verify that the sectors TC± m will not be 
concerned by this operation if \m\ > n. The resulting entanglement will receive contribution 
only from the sectors 7i m with |m| < n and will be a function of A" and n. This can be seen from 
equation (J3.37)) as follows. The first row Oa of the above matrix corresponds to with 
a = m and 0<m<n. The second row (j)\ a corresponds to Q% with a = m + 1 and — l<m<n — 1. 
The nth row n -ia corresponds to with a = m + n — 1 and — (n — l)<m<l while the last 
row na corresponds to with a = m + n and — n<m<0. Thus the unavailable degrees of 
freedom are with —n<m<n and l<a<n + 1. Remark that for any fixed m>0 the index a 
runs over the range \<a<n + 1 — m.For m < the index a runs over the range 1 — m<a<n + 1. 
The resulting reduced density operator is therefore 

n 

Pred = (g) P tl (3-39) 

m=—n 

and the corresponding entropy is given by ( also with = p[ c ^ ) 



S N (n) = S + 2j2Sm (3-40) 

m=l 

where S m is computed using the matrix given by 

n+l— m 

A M = _ £ Wr a ^ W £\ (3.41) 



a=l 



The available indices i, j run from n + 1 — m to n. Again we define a scaled entropy S by 

S = M^. (3.42) 
2n + 1 v ; 

For a fixed we compute S numerically for n = 2 ■ ■ ■ ■, N — 1. The results for A^ = 400 
are depicted on Figure 3. For all values of n with 30 < n < 370 the scaled entropy is almost 
a constant being S = 0.235 for n = 75 and drops very slowly until it reaches S = 0.230 for 
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n = 350. The fact that S can not be a constant is due to the finiteness of ratio of n to N, which 
is related to the fact that the entropy can not be a constant all the way to n = N, since at this 
point we will have traced out all the disc and the result must be zero. When n approaches N i.e 
n > 380 the entropy starts to decrease faster until it reaches zero. Therefore for iV ^> n ^> l 1 , 
we can safely write 

S N (n) = 0.23(2n + 1) (3.43) 
Here again we find the entropy to be directly proportional to the number of degrees of freedom 

0.45 

0.4 
0.35 

0.3 

| 0.25 

if 02 
Oi 

0.15 
0.1 
0.05 


50 100 150 200 250 300 350 400 

n 

Figure 3: The scaled entropy S for a massless scalar resulting from tracing out a sub-disc of a 
fuzzy disc . 



x 



on the boundary, and as can be easily seen there are 2n + 1 of them. It is interesting to note 
that the leading contribution to Sjy(n) is independent of N. 

Let us now consider the fuzzy disc and instead of tracing out a smaller sub-disc consider 
the entropy resulting from tracing out half of the disc (upper or lower half). The procedure is 
exactly similar to the fuzzy sphere case. We define a scaled entropy S = Sn/N and compute 
numerically the resulting entropy for a massless scalar for each N. The results are depicted on 
Figure 4. It is seen that S is rapidly converging towards a constant value independent of N 
equal to 0.341, namely 

S N = 0.3AN. (3.44) 

In view of our previous discussions this result is also expected. The entanglement entropy is 
proportional to the number of boundary degrees of freedom which is iV + 1 in this case. 

lr Th.e condition N 3> n may not be needed to obtain this result, although it is hard to decide by numerical 
calculation how close n can be to N. 
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Figure 4: The scaled entropy S for a massless scalar resulting from tracing out half of a fuzzy 
disc . 

Now, the only difference between equation (|3.43jl and equation (|3.44j) is the coefficient of 
proportionality which can be attributed to the shape of the boundary. When we traced out a 
sub-disc D n the boundary is a sort of a fuzzy circle with 2n + 1 degrees of freedom, whereas 
tracing out the lower or upper half of the disc the resulting boundary is linear with N degrees 
of freedom. 

Before going into the interpretation of our results we note that in the fuzzy disc we have 
restricted the calculation to the massless case. Indeed unlike the fuzzy sphere none of the 
matrices is singular, and this is due to the fact that the fuzzy Laplacian in the massless 

case has no zero modes because of the breaking of translation invariance, therefore the massless 
case does not need special treatment. Furthermore for massive field it can easily be shown 
(numerically) that unless the mass is of the order of the cutoff e.g /i 2 6* ~ 1, the resulting 
entropy has the same value as the massless case. However, if the fi 2 8 ~ 1 the resulting entropy 
differs slightly from the massless case but starts to converge towards the value of the massless 
case for large N. 

Now we turn to the area interpretation of our results. First it must be noted that despite 
the work of [14] the fuzzy disc is still a poor understood object and our discussion are at the 
tentative and qualitative level . The radius of the fuzzy disc is defined through the equation 
R 2 = N9. In order to interpret equation (|3.44|) as an area low in two dimension let us define 

the effective short distance cutoff for the disc 2 to be Ajv = R/N = J j^. Thus equation (|3.44|) 



2 Indeed the entropy is dimensionless and in the free massless theory our action does not know about the 
scale 9 and all we can obtain is a pure number. However the fact that the entropy is proportional to the number 
of degrees of freedom on the boundary makes the area interpretation natural. 
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becomes ( with D = 2R is the diameter of the disc ) 

S N = 0.17^-. (3.45) 
Xn 

In the continuum limit iV — > oo with R fixed ( the limit of the continuum disc of radius R ) we 
recover the usual UV divergence since in that limit — >0. 

Similarly we can write equation ()3.43j) ( with r 2 = n6, \ n = R/y/Nn = \J\ and C is the 
circumference of the disc given by C = 2irr ) as 

S N (n) = — ^. (3.46) 

Again in the commutative limit n — > oo with r fixed ( the limit of the continuum disc of radius 
r ) we recover the usual UV divergence. 

Although the two resulting entanglement entropies considered in the fuzzy disc case diverge 
in the commutative limits of the continuum discs they have different limits when N — > oo with 
9 fixed which is the noncommutative limit of the Moyal plane. 

In the case of equation ()3.46|) if we keep the size of the small disc finite i.e n finite, the 
result remains finite even in the limit N — > oo. Taking our calculation as a regularization of the 
Moyal plane then we would expect the entanglement entropy resulting from ignoring a finite 
region (disc) to be finite and proportional to the number of degrees of freedom on the boundary 
of the ignored region. 

For the case of equation ()3.45|) taking the limit iV — > oo amounts also to recovering the 
Moyal plane. However the region which we are ignoring now is half of the Moyal plane and 
it is not of finite size. From equation (|3.45jl we see that the entropy per unit length which is 

given by S/D = 0.17Wy diverges as N approaches infinity. But this divergence is not of UV 
origin since it shows up as a consequence of blowing up the ignored region, i.e the divergence 
is a consequence of integrating out an infinite number of degrees of freedom. Furthermore we 
observe that regardless of the value of N the entropy per unit length diverges in the limit # — > 
which is the standard UV divergence. Indeed in the continuum case the entropy ( per unit area 
) can be rendered finite by just introducing short distance cutoff in the normal direction, and 
no IR cutoff is needed. Therefore the divergence of the entropy in this Moyal plane limit can 
be understood as a UV-IR mixing rather than just coming from UV origin. 



4 Conclusions and outlook 

We have computed numerically the entanglement entropies of a scalar field on different 
fuzzy spaces. For the fuzzy sphere we have shown that the entropy resulting from tracing out 
half of the sphere is proportional to the number of non-commutative degrees of freedom on the 
boundary. 

The same problem was considered on a fuzzy disc and two ways of tracing were considered. In 
the first case we computed the entropy resulting from tracing out the degrees of freedom residing 
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in a smaller sub-disc. In the second case we considered the entropy resulting from tracing half 
of the fuzzy disc. In both cases the entanglement entropy turned out to be proportional to 
the number of degrees of freedom on the boundary. In the commutative limit both entropies 
suffer from the standard UV divergences. However in the Moyal plane limit the two entropies 
showed different behavior. Whereas the first entropy, resulting from ignoring the degrees of 
freedom inside a fuzzy sub-disc, remained finite when the disc becomes a Moyal plane, the 
second entropy per unite area (length) resulting from ignoring half of the fuzzy disc diverges in 
the Moyal plane limit. 

The divergence in the Moyal limit could easily be seen to be arising from an UV-IR mixing. 
Therefore we would in general expect non-commutativity to render the entanglement entropy 
finite as long as the ignored region is of a finite size, as is the case for the black hole. However, 
if the ignored region is not of a finite size, the entropy per unite area may still be divergent but 
now due to an IR-UV mixing. Finally it would be interesting to investigate the divergences of 
the entanglement entropy per unite area in the Moyal plane using analytical method not just 
numerical. One possible approach would still be to approximate the Moyal plane by a fuzzy 
disc and use 1/N expansion to evaluate the entropy. 
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